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FOREWORD 


This  document  is  Volume  II  of  a  four  volume  report. 
The  four  volume  report  is  titled,  "NAVSTAR  Global 
Positioning  Systems  Special  Studies  and  Engineering 
Program".  Volume  II  summarizes  the  work  performed 
for  the  Air  Force  and  addresses  two  subjects.  The 
first  subject  deals  with  Generalized  Cross-Correla¬ 
tion  Projection  of  GPS  Codes  in  Cases  of  Differential 
Doppler  Shifts ,  and  was  submitted  in  draft  form  to 
SAMSO  on  11  July  1975  as  STI  TR-7115-1. 

The  second  subject  deals  with  the  GPS /TRIDENT  Code 
Design  and  was  submitted  in  draft  form  on  August  6, 
1975  as  STI/GPS-051 . 
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GENERALIZED  CROSS -CORRELATION  PROPERTIES  OF  GPS  CODES 
IN  CASES  OF  DIFFERENTIAL  DOPPLER  SHIFTS 


I.  INTRODUCTION 

In  the  coruse  of  GPS  (Global  Positioning  System)  code  design, 
32  distinct  Gold  Codes  from  the  same  family  of  period  1023 
bits  were  selected  as  the  C/A  (clear /acquisition)  codes,  each 
for  one  of  the  32  NAVSTAR  satellites.  There  are  at  least 
three  main  reasons  for  such  a  selection. 

1.  Gold  code  of  period  1023  bits  is  a  £ >od  ranging 
code:  its  (normalize)  autocorrelation  function 
has  a  maximum  value  (0  dB)  when  compared  with 
itself  and  very  small  values  (■  -23.9, or  -24.2, 
or  -60.2  dB)  when  compared  with  its  m-bit 

(m  f  0)  shifts  such  that  false  code-phase  lock 
may  be  avoided  during  acquisition. 

2.  Gold  codes  from  the  same  family  of  period  1023  bits 
are  good  multiple-access  codes  for  spread-spectrum 
modulation:  they  have  very  small  cross-correlation 
(-  -  23.9,  or  -25.2,  or  -60.2  dB  in  the  case  of 
zero  Doppler  shift)  between  each  other  and  therefore 
can  be  transmitted  simultaneously  in  the  same 
frequency  band  without  severe  mutual  interference. 

3.  The  1023 -bit  period  (-  1  m  sec  at  bit  rate  1.023 
Mbps)  appears  to  be  a  reasonable  compromise  between 
the  correlation  magnitude  and  the  acquisition  time: 
the  code  period  has  to  be  long  in  order  to  achieve 
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small  autocorrelation  (m  ^  0)  and  cross-correlation, 
whereas  quick  acquisition  normally  requires  the  code 
period  to  be  short. 

The  good  cross-correlation  property  of  GPS  C/A  codes  in  the 
case  of  zero  Doppler  shift  has  been  well  known  for  years.  Less 
well  known  are  the  cases  when  there  are  finite  Doppler  shifts 
involved.  Sec.  II  defines  the  problem.  Computer- simulation  re¬ 
sults  on  the  statistics  of  the  interference  levels  are  presented 
in  Sec.  III.  Finally,  in  Appendix  A,  the  generation  of  the  32 
C/A  codes  are  depicted. 

II.  DEFINITION  OF  THE  PROBLEM 

Before  proceeding  with  the  details,  it  is  convenient  to  agree  on 
notation.  Capital  letter  "G,  X,  Y,— -"will  usually  be  used  to 
denote  periodic  sequences,  while  lower-case  letters  "g^,  x^,  y^, 
— M(i  *  1,  2,  3,  — )  will  generally  stand  for  the 
ick  bits  of  the  periodic  sequences.  To  distinguish  between  a 
periodic  sequence  and  its  cyclic  shift,  let,  for  example,  Gj*  G 
be  the  periodic  sequence  starting  with  the  first  bit  of  G: 

-  G  -  gL  g2  g3 -  (1) 

then  G^^  will  designate  the  m-bit  cyclic  shift  of  G,  starting 
with  the  (mfljth  bit  of  G^: 

Gl+m  “  8l+n>  82+m  83+m  *  ‘  ^ 

With  these  preliminaries,  consider  the  simplified  GPS  configur¬ 
ation  in  Figure  1,  where  only  6  of  the  32  NAVSTAR  satellites  are 
shown.  The  C/A  signal  C j ( t)  transmitted  by  the  jth  satellite 
at  time  zero  can  generally  be  expressed  as  a  biphase-modulated 
carrier  given  by: 

Gj(t)  -  Dj(t)  G^kj)  cos  2*fct;  j-lf  2, —  32 

(3) 
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dsaa 


***** 


where 

Dj (t)  =  50  bps  data  bit  stream 
f  =  1575.42  MHz  (L-, )  or  1237.83  MHz  (L9) 

c  1  z  w 


and 


G1(kj)  Xx  ©  Y1+k^ 


A.  L 

is  the  j  C/A  code  generated  by  modulo- 2  sum  of  two  maximal- 

length  sequences  Xi  and  Yi+fc  (see  Appendix  A).  Because  normally 

J 

all  satellites  have  different  distances  and  radial  velocities 
with  respect  to  the  user  receiver,  the  arriving  C/A  code  phases 
(due  to  different  distances  or  time  delays)  and  carrier  frequencies 
(due  to  different  radial  velocities  or  Doppler  shifts)  will  be 
different.  Symbo  lie  ally,  let  C^(t)  denote  the  arriving  C/A 
signal;  then, 


j  -  1,  2,  — ,  32  (5) 

where 


Gm .(kj)  “  (roj-1)  "  cyclic  shift  of  G^(kj) 

Yj  *  the  time  delay  between  the  j ^  satellite  and 
the  user  receiver. 


-  the  Doppler  shift  of  the  jth  carrier 
frequency. 
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and 

m.^l-l.OZS  Mbps  • y .  (modulo  1023) 

J  J 

v  j  =  the  relative  radial  velocity  between  the 
satellite  and  the  user  receiver 
c  =  velocity  of  light 

Because  the  data  bit  stream  usually  cannot  worsen  the  cross¬ 
correlation  properties  to  be  discussed,*  consider  for  simplicity 
of  analysis  only  the  remaining  terms  in  (5).  Then  the  input  to 
the  user  receiver  will  be  the  summation  of  Cj  (t)'s  given  by: 


(6) 


Refer  to  Figure  2,  the  simplified  functional  diagram  of  the  user 
receiver  channel  constructed  to  receive  the  j*"*1  C/A  signal,  and 
suppose,  for  ease  of  discussion,  there  is  only  one  interf erring 
signal  (or  undesired  signal)  C^(t)  present  (which  can  be  readily 
generalized  to  cases  when  there  are  more  than  one  interferring 
signals).  The  reference  signal  Rj(t)  may  be  obtained  from  the 
output  of  VCO  (voltage-controlled  oscillator)  when  the  carrier 
PLL  (phase- locked  loop)  is  locked  onto  C.(t).  The  local  code 

A.L  J 

generator  produces  the  k.  C/A  Code  GR(k.),  and  t*ie  code 

phase  selector  shifts  the  code  phase  n  until  n^mj  (or,  equivalently, 
until  the  power  meter  exceeds  about  0  dB) .  At  this  moment,  the 


t  One  of  the  reasons  is  that  it  has  period  (20  m  sec)  equal  to 
integral  multiple  of  the  period  of  C/A  codes  (1  m  sec)  and  hence 
has  no  effect  on  the  spectra  (or  the  cross- spectra)  of  the  C/A 
codes  to  be  discussed. 
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SIMPLIFIED  FUNCTIONAL  DIAGRAM  OF  THE  RECEIVER  CHANNEL  TO  DETECT 


input  S(t)  to  the  integrator  (or  low- pass  filter)  becomes 


S(t)  =  14 -Gm  (kj)  Gm>  (k^  cos  (2tt  fdt  4-</>) 


where 


£d-Afi-A£j 


is  the  differential  Doppler  shifts  between  the  i*1  and  the  j  C/A 
signals  arriving  at  the  user  receiver.  Note  that  before  the 
desired  code  phase  is  found  (viz.,  before  n  =  iik),  the  5  term 
in  (7)  is  given  by 


Gn(kj)  Gm.(ki)  cos  (2»fdt+$) 
%  (k^)  cos  (2«-fdt  +  <f>) 


where 


W  -  Gn<kj>  Gmi(ki)  ;  k/  *  kj  °r  ki 


is  a  third  Gold  code  from  the  same  family. 


The  contribution  of  8  to  the  power-meter  readout  can  hence  be 


written  as 


P  .  IJ^T  G/(k^)  cos  (2irfdt  +  P)  dt 
“  f  I  Fje  <fd>  I  for  tj  <  1  kHz  and  fd  «.<  1  MHz 
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where 

Fx(f)=?{G/(k £)}  (12) 

is  the  Fourier  spectrum  of  G^(k^).  Since  (^(k^)  is  periodic  with 
period  **  1  msec,  F^(f)  is  a  line  spectrum  having  components  at 
+  N  kHz,  N  being  any  integer. 

III.  CUMULATIVE  PROBABILITY  OF  THE  INTERFERENCE  LEVEL  P 

As  n  in  (9)  shifts  quickly  bit-by-bit  to  search  for  the  code  phase 
°f  Gm.(kj)>  mi  stays  relatively  unchanged;  consequently,  G^(k^) 

becomes  the  whole  family  of  Gold  ''  des.  The  cumulative  probability 
of  the  interference  level  P  can  therefore  be  derived  from  the 
Fourier  spectra  of  the  whole  family  of  Gold  codes,  as  shown  in 
Figure  3-8  for  fd  -  0,  1,  2,  3,  4,  and  5  kHz.  (For  most  GPS  user 
receivers,  the  differential  Doppler  shift  f^  S  5  kHz).  In 
addition,  the  maximum  value  of  P  for  each  case  is  tabulated  in 
Table  2. 


Differential  Doppler 
Shift  between  Two 
Arriving  C/A  Signals 
at  the  User  Receiver 

Max  Value  of  INTERFERENCE 
LEVEL  P 

0  Hz 

-23.9  dB 

1  kHz 

-21.1  dB 

2  kHz 

-21.1  dB 

3  kHz 

-21.6  dB 

4  kHz 

-21.1  dB 

5  kHz 

-21.9  dB 

TABLE  2 

MAXIMUM  VALUE  OF  INTERFERENCE  LEVEL P 
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CUMULATIVE  PROBABILITY 


INTERFERENCE  LEVEL  P(dB) 

FIGURE  5 

CUMULATIVE  PROBABILITY  OF  INTERFERENCE 
LEVEL  FOR  GPS  C/A  CODES 


CUMULATIVE  PROBABILITY 
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INTERFERENCE  LEVEL  P(dB) 

FIGURE  8 

CUMULATIVE  PROBABILITY  OF  INTERFERENCE 
LEVEL  FOR  CPS  C/A  CODES 
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GENERATION  OF  THE  32  GPS  C/A  CODES 


Let  X,  and  Yx  ,  both  of  period 
length  sequences  that  generate 
in  Figure  A.l: 


1023  bits,  be  the  two  maximal- 
the  32  GPS  C/A  codes  as  illustrated 


Gl(kj)  “  xl  ®  Yl+kj 


«.  that  T1+k  *.  V“*  —  Shl£t  °£  Vl>  U  "  ? 
taking  tbe  n»dulo-2  sun  of  and  •,.  .  ,»*  °*  «’*“  '*'} 

generator.  Tbe  values  of  H  -  **  “" 

!,  rh.  resulting  GPS  C/A  nodes  are  showi  In  Table  A. 
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TABLE  A.1 


Code  Phase  Selections  and  the  First  10  Bits 
of  the  Resulting  C/A  Codes 


j 

si 

S2 

First  Ten  Bits  of  [Gn(K^)j 

Decimal 

Binary 

(From  Left  to  Right) 

1 

2 

6 

800 

1100100000 

2 

3 

7 

912 

1110010000 

3 

4 

8 

968 

1111001000 

4 

5 

9 

996 

1111100100 

5 

1 

9 

603 

1001011011 

6 

2 

10 

813 

1100101101 

7 

1 

8 

601 

1001011001 

8 

2 

9 

812 

1100101100 

9 

3 

10 

918 

1110010110 

10 

2 

3 

836 

1101000100 

11 

3 

4 

930 

1110100010 

12 

5 

6. 

1000 

1111101000 

13 

6 

7 

1012 

1111110100 

14 

7  • 

8 

1018 

1111111010 

15 

8 

9 

1021 

1111111101 

16 

9 

10 

1022 

1111111110 

17 

1 

4 

622 

1001101110 

18 

a 

5 

823 

1100110111 

19 

3 

6 

923 

1110011011 

20 

4 

7 

973 

1111001101 

21 

5 

8 

998 

1111100110 

22 

6 

9 

1011 

1111110011 

23 

1 

3 

563 

1000110011 

24 

4 

6 

966 

1111000110 

25 

5 

7 

995 

1111100011 

26 

6 

8 

1009 

1111110001 

27 

7 

9 

1016 

1111111000 

28 

8 

10 

1020 

1111111100 

29 

1 

6 

599 

1001010111 

30 

2 

7 

811 

1100101011 

31 

3 

8 

917 

1110010101 

32 

4 

9 

970 

1111001010 

GPS /TRIDENT  CODE  DESIGN 
1  .  Introduction 

This  report  summarizes  the  results  of  a  GPS/TRIDENT  code  design 
under  contract  to  GPS  (Global  Positioning  System),  Joint  Program 
Office  (F04701- 74-C-0301) .  The  main  purpose  of  this  study  is  to 
devise  a  ranging  code,  called  H  code,  for  trajectory  determina¬ 
tion  of  the  TRIDENT  missile  employing  NAVSTAR  GPS.  In  addition 
to  being  a  good  ranging  code,  the  H  code  is  further  required,  for 
spread- spectrum  multiple-access  (SSMA)  operation  in  GPS/TRIDENT 
system,  to  be  a  good  spectrum- spreading,  multiplexing  code  in 
the  sense  that  it  has  very  small  interference  or  cross-correlation 
(<  -30  dB)  with  any  one  of  the  32  C/A  (clear /acquisition)  codes 
transmitted  by  the  GPS  NAVSTAR  satellites. 

Figure  1  illustrates  the  basic  concept  of  GPS/TRIDENT  svstem. 

All  six  C/A  codes  (which  are  also  good  ranging  and  multiplexing 
codes,  each  from  one  of  the  six  NAVSTAR  satellites)  and  the  H 
code  [from  the  f own  range  support  ship  (DRSS)]  are  transmitted 
synchronously  (vith  atomic  standards  on  all  transmitter  clocks) 
to  the  TRIDENT  fissile  in  the  same  frequency  band  L^.  Because 

all  seven  sources  normally  have  different  distances  and  different 
radial  velocities  with  respect  to  the  TRIDENT,  the  arriving  code 
phases  and  carrier  frequencies  will  be  different  from  one  access 
to  another  and;  furthermore,  from  one  instant  to  another  as  •'ll 
vehicles  are  moving  with  time. 

The  TRIDENT  transponder  translates  the  received  signals  from  Lj- 
to  S-band  by  multiplying  the  received  signals  with  an  auxiliary 
sinusoid  and  then  coherently  retransmits  these  frequency- 
translated  signals  via  the  down- link  to  the  DRSS.  Basically 
each  signal  treats  the  other  six  signals  as  additive  interference 
because  each  signal  will  be  received  by  one  of  the  seven 
independent  channels  of  the  DRSS  receiver.  A  receiver  channel 
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consists  mainly  of  a  phase-locked  carrier  tracking  loop, 
a  delayed-lockeu  code  tracking  loop,  and  a  correlation 
detector.  Theoretically,  the  phase  of  a  desired  code  may  be 
found  by  the  following  steps: 

1.  Cross-correlate  the  received  signals  with 

a  locally  generated  replica  of  the  desired  code. 

2.  Shift  the  phase  of  the  local  code  until  maximum 
correlation,  which  in  the  absence  of  interference 
and  noise  has  a  normalized  value  of  1  or  0  dB 

at  the  correlator  output,  is  reached.  At  this 
moment,  the  local  code  is  said  to  be  locked  onto 
(i.e.,  in  phase  synchronization  with)  the  desired 
code. 

Once  the  phase  and  hence  the  propagation  delay  of  every  received 
code  is  known,  the  position  of  TRIDENT  can  then  be  easily 
calculated. 

Evidently,  to  avoid  false  code-phase  lock  during  acquisition, 
it  is  desired  that  each  of  the  seven  codes  (-  six  C/A  codes 
plus  the  H  code)  has  large  correlation  when  compared  with  it¬ 
self  and  very  small  correlation  when  compared  with  (1)  its 
m-bit  shift  (m  being  any  integer  4  0  modulo  the  period  of 
the  code)  or  (2)  any  of  the  other  six  codes. 

In  the  course  of  GPS  code  design,  32  distinct  Gold  codes  from 
the  same  family  of  period  1023  bits  (-  1  m  sec  at  bit  rate 
1.023  Mbps)  were  selected  as  the  C/A  codes  (see  Appendix  A 
for  their  generations),  each  for  one  of  the  32  NAVSTAR  satellites. 
As  can  be  seen  from  Appendix  B,  all  32  (and  hence  any  subset 
of  six)  C/A  codes  are  good  ranging  md  multiplexing  codes  and 
have  the  desired  properties  stated  in  the  previous  paragraph; 


3 


only  the  H  code  remains  to  be  designed.  Sec.  2  presents 
four  promising  candidates  of  the  H  code  and  their  implementa¬ 
tions.  Based  upon  a  simplified  functional  diagram  of  the 
correlation  receiver,  Sec.  3  analyzes  the  effects  of  inter¬ 
fering  codes  from  other  accesses  at  the  correlator  output. 
Computer- simulation  results  of  the  generalized  correlation 
properties  of  these  four  H  codes  are  briefly  described  in 
Sec.  4;  further  details  of  these  results  may  be  found  in 
Appendices  C  through  E  in  the  form  of  computer  listings.  All 
numbers  in  the  listings  are  in  dB  and  read  from  left  to  right 
then  top  to  bottom.  Section  5  gives  some  general  remarks  on 
Hl/H2-type  codes.  A  CC*  code  (a  code  succeeded  by  its  conjugate) 
which  was  examined  (but  not  recommended)  is  discussed  in 
Section  6. 


2.  FOUR  RECOMMENDED  GPS /TRIDENT  CODES  AND  THEIR 

IMPLEMfeNTATIOhS  — 


Before  proceeding  with  the  details,  it  is  convenient  to  agree 
on  notation.  Capital  letter  "G,  H,  X,  Y.  — "  will  usually 
be  used  to  denote  periodic  sequences,  while  lower-case  letters 
"®i*  ^i*  xi»  Yi>  ' — "(i  “  1*  2,  3,  — )  will  generally 

stand  for  the  iC^  bits  of  the  periodic  sequences.  To  dis- 
tinquish  between  a  periodic  sequence  and  its  cyclic  shift, 
let,  for  example,  G^  »  G  be  the  periodic  sequence  starting 
with  the  first  bit  of  G: 

Gx  “  G  “  g2  g3  -  -  -  (1) 

then  Gj^  will  designate  the  m-bit  cyclic  shift  of  G,  starting 
with  the  (nri-l)*  bit  of  G^ 

Gl+m  -  ®l+m  ®2+m  «3+m -  <2> 

Further,  unless  otherwise  specified,  the  following  conventions 
will  be  used  throughout  this  report: 
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X,  Y  =  the  two  psuedo-noise  (PN)  sequences  that 
generate  the  32  GPS  C/A  codes 

G  (kj)  -  the  jth  GPS  C/A  code  (j  =  1,  2,  — ,  32) 

H  -  the  GPS /TRIDENT  code 

and  for  convenience  of  discussion,  the  four  recommended 
GPS/TRIDENT  codes  will  be  denoted  by  HI,  H2,  H3,  and  H4. 

Figure  2-5  depict  the  implementations  of  the  four  GPS/TRIDENT 
codes  with  their  parameters  of  interest  shown  in  Table  1. 


Explicity, 

H1  -  X1  X2  X3  '  -  -  X32 

(3) 

H2  “  X1  X2  X3  •  -  *  X64 

(4) 

are  sequences  of  cyclic  shifts  of  X ,  and  H3  and  H4  are  short¬ 
ened  (or  short-cyled)  PN  sequences: 

H3  -  a  periodic  sequence  having  period  equal  to  the 
first  32  x  1023  bits  of  a  PN  sequence  of  period 

215  -1  (-  32  x  1023  +31)  bits 

H3  ■  a  periodic  sequence  having  period  equal  to  the 
first  64  x  1023  bits  of  a  PN  sequence  of  period 

216  -1  (-  64  x  1023  +63)  bits 

To  further  manifest  the  short-cycling  operation,  Figures  6 
and  7  illustrate  the  state  diagrams  of  the  H3  and  H4  generators; 
the  initial  state  (or  Sate  1)  is  defined  as  the  state  when 
the  contents  of  the  sequence  generator  are  "all  ones11. 


Because  both  the  32  and  64  m  sec  clock  pulses  (which  control 
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ED  GPS /TRIDENT  CODES 


the  periods  or  short-cyclings  of  the  «c^esb^ 

the  rr/sny  ^rfr”3  Mbps^chip  rate  pulses 

riri'Ces ». 

frequency-standard  t0  be  equal  to  integral 

four  recommended  H  codes  the  GPS  C/A  codes, 

multiples  of  the  period  (  1023  '  the  subsequent 

thereby,  to  V^rost-cowelation  between  the  H  codes 
Tthe  GPS  CM  codes.  Finally,  it  is  worth  mentioning  that 
If  the  32  and  64  n  sec  cloch  pulses  in  Figures  4  and  are 
removed ,  then  the  H3  and  «4  of  feedback 

r;i— 

215  -1  and  21&  -1  bits,  respectively. 


3.0 


C ation  Receiver  Analysis 


According  to  the  discussions  presented  in  Section  1.  the  H  code 

ference^  consider  ^Figure  7TslTu?iei  functional  di  agram of 
the  DRSS  receiver  intended  to  detect  the  H  code.  Let  theint 
faring  signal  C.(t)  b.  only  on.  of  th.  .1.  CM  sign* •  and  o  th. 

youar  ..  th.  daslrad  signal  D(t).  If  for  .Inyllotty  of 
analysis  th.  absolute  tin.  and  coda-phas.  tafotsocas  ar.  naglaat 
(which  will  not  change  the  result  we  are  going  to  derive) .  then 
the  input  to  the  receiver  can  be  simply  written  as 


S(t)  -  D(t)  +  Ca(t) 

-  Hj.  cos  2v  fHt  +  G1+ffi  (kj)  cos  2n  (fH  +  tjt 


(5) 
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SIMPLIFIED  FUNCTIONAL  DIAGRAM  OF  THE  RECEIVER  CHANNEL  TO  DETECT 

H  CODE 


where 


fH  *  carrier  frequency  of  D(t) 

m  »  the  relative  code  phase  between  the  H  code  carried  by 
D(t)  and  the  jth  C/A  code  G  (kj)  carried  by  Cj(t) 
fd  -  the  differential  Doppler  shift  or  the  carrier-frequency 
difference  between  D(t)  and  Cj (t) 

Note  that  m  and  f^  are  caused  by  the  fact  that  both  sources  have 
different  distances  and  radial  velocities  with  respect  to  the 
TRIDENT,  and  not  by  coherent  retransmission  from  the  TRIDENT  to 
the  DRSS.  As  sources  are  moving  with  time,  m  and  f^  will  also 
change  from  one  instant  to  another. 

When  both  the  phase-locked  loop  (PLL)  and  the  local  H  code 
generator  are  locked  onto  D(t),  the  input  to  the  finite-time 
integrator  becomes 

S(t)  -  1  +  Hx  G1+  m  (kp  cos  2  it  fdt  (6) 

Because  the(0  to  T)  finite -time  integrator  can  normally  be 
approximated  as  a  low-pass  filter  of  bandwidth  ip,  the  contribution 
of  the  interference  term  6  in  (6)  to  the  power-meter  readout  can 
hence  be  expressed  as 

|FGH<kj*  m*  fd>|  (7) 

where 

jFGH  (kj*  f>  -  |  ^|H1  G1  +  m 

is  defined  as  the  cross-spectrum  between  and  +  m  (kj). 
m  being  the  relative  code  phase  between  H  and  G  (kj)* 
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As  evident  from  (7),  as  m  and  are  changing  with  time,  the 

requirement  of  small  interference  level  P  at  the  correlator 

• 

output  is  equivalent  to  the  requirement  of  small  cross-spectrum 
between  and  for  all  possible  relative  code  phases 

m  =  0,  1,  2,  -  -  -  1022  bits 

and  all  possible  differential  Doppler  offset  f^  (ranging  from  0 
to  5  kHz  in  most  cases  of  practical  interest) .  (As  will  be  seen 
later  in  Sections  4,2  and  4.3,  all  four  recommended  H  codes  appear 
to  haveP<-.30  dB  required.) 

4.0  GENERALIZED  CORRELATION  PROPERTIES  OF  GPS /TRIDENT 

CODES 

4.1  Autocorrelation  Properties 

Define  the  autocorrelation  function  of  the  H  code  as 


A  (m)  -  -  E  h4  h. . 

P  i_i  i  1+in 


where 


the  period  of  the  H  code ,  bits 


h^  -  the  icn  bit  of  the  H  code  *  1  or  -1 

hi+m  -  the  (i  +  m)th  bit  of  the  H  code  -  1  or  -1 

m  *  the  relative  code  phase, bits 


then  it  is  desired  that  (Section  1) 


A(m)«  A(0) 


for  any  m  +  0  (module  p) .  All  four  recommended  H  codes  have  such 
a  desired  property  as  evidenced  in  Figures  9-12,  where  A(m)  in  dB 
is  defined  as  20  log^Q  A(m)  and  only  cases  when  m  -  0  to  1000  bits 
are  shown  because  the  receiver  normally  will  have  some  knowledge 
about  the  position  of  the  TRIDENT,  say  within  accuracy  of  +300  km 
(«Hl000-bit  code-phase  delay  at  chip  rate  1.023  Mbps)  or  so. 


Figure  9  Autocorrelation  Function  of  H  1 


Figure  10  Autocorrelation  Function  of  H  2 


0  dB 
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dB 


Detailed  listing  of  these  four  autocorrelation  functions  may  be 
found  in  Appendix  A. 

4#2  Cross-Correlation  Properties  in  the  Case  of  Zero 

Differential  Doppler  Offset 

In  the  special  case  when  the  differential  Doppler  offset  fd  =  0 
the  interference  level  P  in  (7)  becomes  the  cross-correlation 
function  Cq^(iii)  be  tween  and  WV 


CGH^  p  6i+m 


p  =  1023  xn  v 

T  =  n  x  1  msec  ' 

n,  ibeing  any  positive  integers.  Computer-simulation  study 
indicates  that  for  H  =  HI 

max  (m)J  S  -33.8  dB  * 

* 

for  all  m  -  0,  1,  . ,  1022  and  all  32  GPS  C/A  codes,  and  for 

H  -  H2  the  corresponding  upper  bound  is  given  by 


CouoWl  *  -37.3  dB 


Further  details  of  these  results  are  shown  in  Appendix  D. 
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4.3  Cross -Spectra  Between  GPS/TRIDENT  Codes  and  GPS  C/A  Codes 

As  was  demonstrated  in  Section  3,  the  interference  level  P 

for  f ^  f  0  is  equal  to  the  cross-spectrum  ^  |FGH^kj’  m’  f^l 

at  f  *  Although  FFT  (fast  Fourier  transformation)  subroutine 

had  been  employed  to  reduce  the  computation  time  required  during  the 

simulation  study  of  the  cross-spectrum,  it  still  takes  about  0.35  min 

CPU  time  of  IBM  360/67,  for  each  m  and  j  selected.  For  each  H  code 

there  are  a  total  of  1023  m's  (m  -  0,  1,  . ,  1022)  and  32  j's 

(j  *  1*2,- — ,32),  or  a  total  of  1023  x  32  cross-spectra  to  be 
calculated,  thereby  requiring  about 

0.35  x  1023  x  32  a  104  min 

CPU  time  to  run  all  possible  cases.  Because  of  such  a  time- 
consuming  nature,  only  several  randomly  chosen  cases  are  studied 
with  the  results  summarized  in  Table  2.  All  four  recommended  H  codes 
appear  to  have  maximum  magnitudes  of  cross-spectra  (or  maximum 
interference  level  p)  less  than  about  -35  dB.  Details  of  these 
results,  again,  may  be  found  in  Appendix  E.  Note  that  all  cross - 
spectra  are  line- spectra  with  frequency  difference'  between  adjacent 
lines  equal  to  31.25  Hz  for  H  *  HI  and  H2  or  15.625  Hz  for  H  -  H2 
and  H4.  To  give  a  better  feeling  about  the  characteristics  of 
these  spectra,  two  particular  spectra  are  illustrated  in  Figures 
13  and  14. 
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H  Codes 


Relative  Code  Phase 
m,  bits 


Maximum  Magnitudes 
of  Cross-Spectra,  dB 


5, 


Some  Remarks  on  Hl/H2-Type  Codes 


An  interesting  property  of  the  Hl/H2-type  codes  is  that  the 
maximum  value  of  the  interference  level  P  decreases  as  the  code 
period  increases.  Stated  another  way,  the  Hl/H2-type  codes  have 
very  good  partial  cross-correlation  property.  To  demonstrate 
this  point,  Figure  15  plots  Pmax  the  maximum  value  of  p  (including 
cases  of  any  differential  Doppler  shifts  less  than  ~  100  kHz)  in 
dB  versus  the  period  p  of  the  codes  from  1  x  1023  to  512  x  1023  bits 
Again,  because  of  the  excessive  amount  of  computation  time  required 
for  an  exhaustive  computer- simulation  study,  this  curve  is  derived 
from  only  limited  computer  runs.  Based  upon  the  statistics  of  these 
limited  computer  runs,  it  is  believed,  however,  that  evan  an 
exhaustive  study  would  not  yield  P  more  than  ~  3  dB  higher  than 
what  is  predicted  by  Figure  15.  A  -3.5  dB/octave  dotted  line  is 
also  shown  for  comparison.  As  can  be  seen  clearly  from  this 
figure,  for  code  periods  between  4  x  1023  and  512  x  1023  bits, 
every  increase  of  the  code  period  by  a  factor  of  2  decreases  Pmflv 
by  a  factor  of  ^3.5  dB 

Another  important  property  of  the  Hl/H2-type  code  which  should 
be  pointed  out  is  that  its  autocorrelation  functions  A(m)  has 
spikes  at  code  phase  shifts  m  -  n  x  1022  and  n  x  1022  +  p/1023 
bits,  where  n  ■  1,2,3,--,  and  p  is  the  period  of  the  code.  As 
an  example,  Table  3  lists  the  spikes  of  A(m)  of  H2.  Notice  that 
A(m)  is  symmetrical  with  respect  to  m  ■  |  32  x  1023  for  H2). 

To  avoid  ambiguities  created  by  these  spikes  during  code  phase 
tracking,  it  is  required  that  the  receiver  has  prior  knowledge 
of  the  received  code  phase  within  accuracy  of  +  1022  bits.  Note 
that  a  similar  (but  not  the  same)  ambiguity  also  exists  in  the 
autocorrelation  function  of  a  single  Gold  code  of  period  1023  bits. 


22 


Table  3  SPIKES  IN  THE  AUTOCORRELATION  FUNCTION  A(m)  of  H2 


Relative  Code  Phase 
m,  bits 

Autocorrelation  Magnitude 
|A(m )|,  dB 

1  x  1022 

-0.1 

1  x  1022  +  64 

-32.3 

2  x  1022 

-0.3 

2  x  1022  +  64 

-28.3 

3  x  1022 

-0.4 

3  x  1022  +  64 

• 

-25.5 

• 

• 

31  x  1022 

• 

-6.0 

31  x  1022  +  64 

-6.5 

32  x  1022 

-6.3 

32  x  1022  +  64 

-6.3 

33  x  1022 

-6.5 

33  x  1022  +  64 

• 

-6.0 

• 

• 

• 

61  x  1022 

• 

• 

-25  5 

61  x  1022  +  64 

-0.4 

62  x  1022 

-2 1.3 

62  x  1022  +  64 

-f  .3 

63  x  1022 

-J2.3 

63  x  1022  +  64 

-0.1 

6. 


A  CC*  Code 


For  a  code  C  ■  —  c  consisting  of  p  bits,  define 

^  ic  it  yc  k 

its  conjugate  c  *  ci  c2  c3  """cp  as  fo^ows: 

1  if  ct  -  -1 
-1  if  c±  -  1 

"k  *-u  ^ 

where  and  are  the  i  bits  of  C  and  C.  respectively. 
Further,  let 

*  *  *  *  * 

CC  *  c^  C2  c3— cp  c^  C2  cp 

be  a  periodic  sequence  of  period  2p.  Then  it  can  be  easily  proved 
k 

that  CC  always  has  zero  (or  -  »  dB)  cross-correlation  with  any 
periodic  sequence  of  period  p. 

In  the  case  of  no  differential  Doppler  shift,  the  interference 

P  between  CC*  and  G  (a  GPS  C/A  code) ,  at  the  output  of  the 

coherent  correlation  receiver,  is  just  the  cross-correlation 

k 

(or  the  DC  term  of  the  cross -spectrum)  between  CC  and  G,  which 
is  -*  dB  for  any  C  having  period  p  -  1023  bits.  Therefore, 

k 

in  this  particular  case,  a  CC  -type  code  would  be  an  ideal 

choice  for  the  GPS/TRIDENT  code.  However,  in  practical  cases, 

there  usually  exists  certain  differential  Doppler  shifts  among 

various  SSMA  channels  such  that  the  magnitude  of  P  becomes  the 

magnitude  of  the  spectral  line  of  the  cross-spectrum  between 

CC  and  G,  at  a  frequency  nearest  to  the  differential  Doppler 

shift.  Unfortunately,  as  one  would  expect  from  the  structure 
k  k 

of  the  CC  code,  the  cross -spectrum  between  CC  and  G,  denoted 

by  F(f),  has  relatively  spiky  magnitudes  at  odd-numbered  lines, 
i.e.,  at  frequencies  f  -  —  Hz.  (Note  that  F(f) 

is  a  line  spectrum  with  frequency  difference  between  adjacent  lines 
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equal  to  ^  -  0.5  kHz).  For  example,  when  C  -  G  -  G^)  and  the 
relative  code  phase  between  CC  and  G  is  m  -  1  blJ* 
simulation  results  indicate  that  the  magnitudes  of  the 
few  spectral  lines  of  F,  in  dB,  appear  as  follows: 


-  oo ,  -28.5,  -78.4,  -22  7,  -72.5,  -31.2,  -68^.6,  -33.1  dB, 

olz  o!skHz  l  1.5kHz  2  kHz  2.5kHz  3kHz  3.5kHz 


As  a  consequence,  CC*  is  not  recommended  to  be  a  good  candidate 
of  the  GPS/TRIDENT  code. 
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